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We study chaotic inflation in the context of modifled gravitational theories. Our analysis covers 
models based on (i) a field coupling ijj{(t>) with the kinetic energy X = — (1/2) g'^'' 8^4)8^4' and 
a nonmimimal coupling (^(j>^R/2 with a Ricci scalar R, (ii) Brans-Dicke (BD) theories, (iii) Gauss- 
Bonnet (GB) gravity, and (iv) gravity with a Galileon correction. Dilatonic coupling with the kinetic 
energy and/or negative nonminimal coupling are shown to lead to compatibility with observations 
of the Cosmic Microwave Background (CMB) temperature anisotropies for the self-coupling inflaton 
potential V{(l>) — A(^*/4. BD theory with a quadratic inflaton potential, which covers Starobinsky's 
! f{R) model f{R) = R-\-R^/{6M^) with the BD parameter ij^BD = 0, gives rise to a smaller tensor-to- 

scalar ratio for decreasing ujbd- In the presence of a GB term coupled to the fleld (f), we express the 
scalar/tensor spectral indices ris and nt as well as the tensor-to-scalar ratio r in terms of two slow- 
, roll parameters and place bounds on the strength of the GB coupling from the joint data analysis 

OA ' of WMAP 7yr combined with other observations. We also study the Galileon-like self-interaction 

$((/))Xn(/) with exponential coupling oc e*"*. Using a CMB likelihood analysis we put bounds 

on the strength of the Galileon coupling and show that the self coupling potential can in fact be 
made compatible with observations in the presence of the exponential coupling with > 0. 
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I. INTRODUCTION 



o 
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2 ■ Inflation, which was originally proposed by a number of authors independently in the early 1980s [J 0|j is at 
, present the main theoretical framework employed in describing the early universe evolution and accounting for the 
■ observational data, specially the observed spectrum of primordial perturbations 0]. The sustained success of this 
framework over the last 3 decades has been impressive, particularly given the enormous improvement in the accuracy 
and resolution of the cosmological data over this period The simplest and most common models of inflation 

considered so far have employed a single scalar field, minimally coupled to the curvature and possessing a canonical 
kinetic term (see @ for reviews). As a result, until recently, much effort has gone into the study of such models, an 
QO ' important example of which has been the chaotic inflationary model [l^l ■ 

, Despite its successes in accounting for important features of observations, however, there is no unique mechanism 
\l • which underpins inflation. Indeed almost since its inception it has been known that an accelerated phase of cosmic 
evolution could be produced by a wide range of mechanisms such as f{R) theories (see the reviews [ill [Si and 
. references therein). Thus, an important task in cosmology has been to narrow down the range of possible alternatives 
and ultimately to situate inflationary models within fundamental theories of physical interactions. There have been 
two approaches to this problem. The first aims to construct individual models that arc directly suggested by super- 
symmetric theories. For ex amp le, chaotic inflationary models have been constructed in the framework of supergravity 
[iSj or superstring theory [ij]. The second, on the other hand, considers classes of generalized models of inflation 
which possess ingredients motivated by field theories such as string theory [l5l - l23 |. 
[ In the absence of a unique, fully successful and non-fine tuned model of the first kind so far, a great deal of effort 
" ■ has recently gone into the study of the models of the second kind. In general such models are expected to possess a 
number of ingredients motivated by fundamental theories, including (a) nonminimal couplings of the field to the Ricci 
scalar R, (b) non-canonical kinetic terms, and (c) higher derivative quantum gravity corrections in their actions, such 
as the Gauss-Bonnet term 

g = R^- AR^pR°^P + Ro^p^sR"^^^ , (1) 

where R is the Ricci scalar, R^p is the Ricci tensor and Rap-ys is the Riemann tensor, or a nonlinear field-interaction, 
for example in the form 

G(0,x)n^, (2) 

where G is in general a differential function of the field (f> and X = —{l/2)g'^^d^(j)d,y4'. In addition, such models would 
also be expected to possess multiple scalar fields, but in order to separate the effects of different ingredients we shall 
confine ourselves to models possessing a single scalar field. 
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In standard chaotic inflation the self couphng A for the inflaton potential V{4>) = A0^/4 is constrained to be small 
(A « 10^^'^) from the WMAP normalization Moreover, this model is in tension with the observations of CMB 
temperature anisotropics because it predicts a large tensor-to-scalar ratio (r w 0.3) 0- If we take into account the 
nonminimal coupling it is possible to realise larger A compatible with the natural values appearing in particle 

physics (A = 0.01-0.1). In the limit where the negative nonminimal coupling ( satisfies the condition \(\ ^ 1, the 
tensor-to-scalar ratio r reduces to the order of 10^'^ with the scalar spectral index rig « 0.96 [l^. This agrees well with 
the CMB observations [26l - [28| . Recently, there has been renewed interest in nonminimally c oup led inflation models 
by identifying the infiaton as a Higgs boson appearing in standard model of particle physics [l7l . [29| . 

It is of interest to see whether the chaotic inflationary models that are in tension with observations can be rescued by 
taking into account the various field couplings mentioned above. In the presence of the nonlinear kinetic interaction 
(l/M^)Xn0 that respects the Galilean symmetry — > -I- 6^ in Minkowski spacetime, for example, it was 
recently shown that the tensor-to-scalar ratio can reduce to r ~ 0.18 for the potential V{(j)) = A0''/4 [sO]. Moreover, 
the coupling constant A is related to the mass scale M by the WMAP normalization, which allows for the possibility 
of having a natural coupling of the order of A = 0.01-0.1. Of course this outcome would be expected to be different if 
we choose more general functions of G{4),X) that depend on both and X. 

In this paper we wish to make a detailed and unified study of infiation in the context of modified gravitational 
theories. To this end, and for concreteness and comparison with previous works, we shall study chaotic infiation, 
sourced by potentials of the form 



where Vq and p are real constants and ~ 1/\/8ttGn = 2.44 x 10^* GeV is the reduced Planck mass (^G N is 
the gravitational constant). We shall consider a number of field couplings such as (i) the non-canonical kinetic term 
uj{(l))X as well as the nonmimimal coupling C,(fpR/2, (ii) Brans-Dicke (BD) theories having explicit couplings (pR and 
{ujbd/4')^ (including f{R) gravity), (iii) Gauss-Bonnet (GB) coupling of the form ^qc'^'^^^p^Q, and (iv) the generalized 
Galileon couphng of the form (e'^"^/*^pVM4"-i)X"n(/) (which reproduces the pure Galileon term [22|, [3l| in the limit 
that /i — > and ?i — > 1). The terms of the forms (i), (iii), and (iv) appear as a next order correction to the tree- level 
action in low energy effective string theory (32j . 

In each model we evaluate the three infiationary observables: (a) the scalar spectral index rig, (b) the tensor spectral 
index nt, and (c) the tensor-to-scalar ratio r. We place observational constraints on the model parameters by carrying 
out a CMB likelihood analysis. We find that in most cases it is possible for the chaotic infiationary potentials with 
p = 2 and p = 4 to be consistent with the current observations. We also show that the equilateral nonlinear parameter 
/nl"'' describing the scalar non-Gaussianity is smaller than the order of unity, apart from in the generalized Galileon 
model with n ^ 1. 

The structure of the paper is as follows. In Sec. |ll]we derive the background equations of motion for the general 
action (|4]) and introduce a number of slow-roll parameters. In Sec. lIIII we present the power spectra of scalar and tensor 
perturbations derived under the framework of linear cosmological perturbation theory. The formula of the equilateral 
non-Gaussianity parameter /^l"'' is also given there. In Sec. llVl we show the three infiationary observables Us, nt, and 
r in the Einstein frame for the theories without the GB or Galileon terms, which are convenient for the analysis in 
Sees. Ivl and IVII Sec.|V]is devoted to the study of the nonmimimal coupling C,4>^R/2 as well as the coupling e^'^/*^p'X. 
In Sec. IVII we consider BD theories with the two potentials V = Vb((/)/Afpi)^ and V{4)) = Vb(0 — Mpi)^ in the Jordan 
frame for arbitrary BD parameters wbd- For p = 2 and wbd = the latter potential covers the Starobinsky's f{R) 
model f{R) = R + R'^ /{6M^). In Sec. I VIII we study chaotic inflation in the presence of the exponential GB coupling 
and place observational constraints on the strength of the GB coupling. In Sec. I Villi we show that the exponential 
Galileon coupling can lead to the consistency of self-coupling chaotic inflation with the observational data. Sec. llXl is 
devoted to our conclusions. 



(P > 0) , 



(3) 



II. MODELS AND BACKGROUND EQUATIONS 



We start with the generalized action 




g -^F{(P)R + cj{cf>)X ~ 1/(0) - a^)g - G(</), X)n0 



(4) 



where 5 is a determinant of the space-time metric 5^^, and is a scalar fleld with a kinetic term X. The functions 
F{(j)), uj{4>), and ^{(f)) are differentiable functions of 4>, whereas G{(j), X) depend on both (f> and X. The fleld (f) couples 
to both the Ricci scalar R as well as the Gauss-Bonnet term Q. 
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We consider the flat Fricdmann-Lemaitre- Robertson- Walker (FLRW) space-time with a scale factor a{t), where t 
is cosmic time. The background equations are then given by 

El = UllFH^ + MllHF -uX -V - 2AH^i - GH^XG.x + 2XG,0 = , (5) 



E2 = MllFH"^ + 2Ml^HF + 2MlyFH + MlyF + ujX-V - IGH^i - 16HHC - 8H^^ - G^x'pX - G.^^^ = 0,(6) 



S3 = {uj + &H^G.x+&Hd^XG^xx-2XG^4,x-2G^4,)'(i) 

+ {3ujH + + 9H^^G^x + SH<j)G^x + SH<P^G^^x - QHG^^ - G 



uj,^X + V^~ &MlH^F^^ - U'llHF^^ + 24.hX^ + 2AH^Hi^ = , (7) 



where H = a/a is the Hubble parameter, a dot represents a derivative with respect to t, and a comma represents a 
partial derivative in terms of cj) or X. Only two of the above equations are independent due to the Bianchi identities, 
ijyEs + Ei+ 3H{Ei - E2) = 0. The combined equation, {E2 - Ei) / {M^^H'^ F) = 0, gives 

_ H _ F F LuX AHi 8Hi 4^ 3^XG,x 4>XG^x 2XG,^ 

^ " "iP - ^2HF ^ 2H^F ^ AI^^H^F ^ M^^F ^ M^^HF ^ M'^^F ^ M^^HF ~ M^^mF ^ AI^^mF ' ^' 

Since e ^ 1 during inflation, the modulus of each term on the r.h.s. of Eq. ([8]) is much smaller than unity (unless 
some cancellation occurs between those terms) . We introduce the following slow-roll parameters 

, F ^ _ iuX .Hi ^ _ ^XG,x , _ 4> , _ XG,4, 



by which we have 



From Eq. ([8]) we obtain 



SpiSp+VF-e), -jj^f; = Sd^F + Vi + ^) ■ (10) 



H^F " ' ^' ' M^,F 



25 X -Sf + 8S^ + 65gx - 4(5g0 + SpiSp + Vf)- SS^Sf + 2S4,Sgx 

2 + 5f-85^ ^ ' 

= Sx -SF/2 + 46^ + 3SGX-25G4. + Oie^), (12) 
where in the latter step we have taken the leading-order contribution. 



III. COSMOLOGICAL PERTURBATIONS 



Let us consider cosmological perturbations about the flat FLRW background. We take into account, up to a 
gauge choice, both the perturbations in the scalar field S(f) and in the scalar and tensor modes of the metric. For 
the calculations of observables, including primordial non-Gaussianities, it is convenient to employ the 4-dimensional 
ADM perturbed metric [s^ of the form 



[(1 + af - a-^{t)e-^'^{dipf] + 2d.,i}dtdx^ + a^{t) (e^^^y + h^j) dx'dx^ 



(13) 



a) and the shift vector di'tp, and it^j 



I] 5 



hi 

to fix the 



where TZ is the curvature perturbation, a and ■0 are related with the lapse (1 

are tensor perturbations. In the metric (|13p we have gauged away a field E appearing as E^ij inside h, 
spatial components of a gauge-transformation vector We choose the uniform- field gauge where 6(f) = 0, in order 
to fix the time-component of ^'^ H^l . 

Expanding the action (|4]) up to second order for the metric (jl3p . performing integration by parts and using the 
Hamiltonian and momentum constraints to eliminate the contribution coming from a and V', we obtain the following 
second-order action 133 



S2 = 



/ 



dt d^xa^C 



7^2 _ 1| (57^)2 



(14) 
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where 



_ wi{AwiW3 + 9wl) 



^2 

3{2wlw2H — W2W4 + 4wi'WiW2 — 2wfw2) 

Wi{4:Wi'W3 + 9w|) 



(16) 



and 



wi EE M'^^F-SH^, (17) 

W2 = M^i{2HF + F) ~24:XG^x -2AH^i, (18) 
W3 = -9M^^FH^ ~9M'^^HF + 3ujX + UAH^i+l8H^(2XG,x+X^G,xx)~Q{XG,^ + X^G,^x), (19) 

W4, = M^iF~8i. (20) 

In order to avoid the appearance of ghosts and Laplacian instabihties we require that 

Q>0, c2>0, (21) 
respectively. One can express (i = 1, • • • , 4) in terms of the slow- roll parameters. For example one has 



W3 = ~9Ml^FH^ (i + 6f- ^5x - l&h - ^^gx + |<5g0 - 25gx\gx + ^Sg^^g^ ) , (22) 



where 



, _ XG^xx , _ XG^^x .„„s 
Aga" = , Ag0 = —p; ■ {26} 

The quantities Xqx and Ag0 arc not necessarily small. 
The expansion in terms of the slow-roll parameters gives 

^2 (5,Y + 4(5gx - 2i5g0 + 2(5g0Ag0 ^^^^ 

Sx + 6Sgx - 2(5g0 + 6(5gxAga: - 25g4,\g<p ' 

= T75^ = + 4(5gx - 2(5g0 + 25g^\g<p 

~2Sg^Sf>^g<p + '^6Sg4,S^\g4, + 26^6gx\gx + 4(5g0'^ga'Ag0 + SSp/A - 12S^Sf + 2SgxS^ - SpSx 
-5SfS Gx + 25FSG<ji + 85^5x + 40(5^(5gx — 16,55(5g0 - 45gai5g,/, + 25gx5x + 48(5| + IS^x 
+0(e3), (25) 



where in the expression for we have picked up the leading-order contributions. In standard slow-roll inflation with 



F = 1, cj = 1, ^ = 0, and G = we obtain the exact expressions = 1 and QjM'K = 8x ~ £• Equation ([M]) shows 



that the nonmininial coupling F{(j))R and the Gauss-Bonnet term (,{4))G do not give rise to contributions to at 
linear order. The effects of those terms on appear at the next order. 
The power spectrum of the curvature perturbation is given by (ssj 



-p ^ ^ li f26) 



which gives the scalar spectral index 



din fc 



-2e - (5q - 3s (27) 
-2e-SF-Vs-s, (28) 



where 

^ _ Q _ Cs _ e 



HQ' ^-Hcs' '"~He/ ^^^^ 
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Wc have assumed that both H and Cg vary slowly, such that dlnfc at c^fc = aH may be approximated by dlnfc 
din a = Hdt. 

The tensor power spectrum is given by [ssj 



Vt 



27r2Q,c? ' 



(30) 



where Qt — wi/A ~ M^^F{1 — 8(5^)/4 and Cj = wa/wi = 1 + 8(5^ + 0{e^). Taking the leading-order contribution in 
"Pt, it follows that Vt c:^ 2H'^ / {tt'^ M'^^F) . The tensor spectral index is 



dink 



(31) 



Ctk—aH 



which is valid at first order in slow-roll. At times before the end of inflation (e <C 1) when both Vg and Vt remain 
approximately constants, we can estimate the tensor-to-scalar ratio, as 



Vs 



16 



Q± 



16ccec . 



(32) 



The non-Gaussianities of scalar perturbations for the action (U) have been evaluated in Ref. [Sg (see also 
Rcfs. [33.l36l.l37j for related works). Under the slow-roll approximation the nonlinear parameter /^l"'' the equilateral 
configuration is 



fequil ^ ( 1 

JNL — 00^ I ^ 



85 
324 
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^ 36 c2 


12 c? 


54 c? 
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10 , 










(-1) 


+ Sgx 



20(1 + Agx) 
81e, 



65 



162c2e^ 



(33) 



where 



A = f2 

s 



(j}H{XG,x + ^X'G,xx + 2X''G,xxx) - 2{2X'G,^x + X'G,^xx)/'S\ , 



12Mp4 



(34) 
(35) 



In the absence of the Galileon term {6gx = = (^00) one has — 1 and Es ~ 5x from Eqs. ((24)) and ([25]) at linear 
order in slow-roll. In this case, the expansion of up to second order gives 



2S^iSF-SS^)i36F-24,5^-ASx) 
Ox 



(36) 



which shows that the GB contribution can only lead to small changes to the value cj. = 1. Then the nonlinear 
parameter in Eq. p3p is approximately given by 



55 
36^ 



5_ 

12' 



10 , 



(37) 



which means that the non-Gaussianity is small for the theories with G = 0. However, the presence of the Galileon 
term can potentially give rise to large non-Gaussianities. 

IV. THE ACTION IN THE EINSTEIN FRAME 

We start by considering non-minimally coupled theories in the absence of the GB and Galileon terms (^ = = G), 
i.e. with actions of the form 



5* 



d x\J —g 



Ml 



F{4>)R + Lo{(f>)X -V{^) 



(38) 
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Since = 1 and s = in these theories, it follows that 



Us - I 



-2e - (5(3 = -2e - (5^ - ?7s — -2es - rjs , 

Q 



16 



16e. 



-8nt , 



where 



{2HF + Ff 



(39) 
(40) 

(41) 



(42) 



In the last approximate equalities of Eqs. ((39|) . pO)) . and (|4T|) we have used the relation eg — e + (5f/2 valid at linear 
order in slow roll. This follows from Eqs. and i.e. e ~ (5^ — i5f/2 and ~ respectively. 

It is convenient to transform the action psp . expressed in the so called Jordan frame, into the one having a scalar 
field minimally coupled to gravity (the Einstein frame), via the conformal transformation 



The transformed action is given by [3S 



Se 



where a hat represents the quantities in the Einstein frame, and 

V 



u = 



p2 ' 



x = / i?(</>)d0, 



'3 ^AfpiF, 



F 



The following relations hold between the variables in the two frames: 



dt 



Defining the variables 



we obtain |29, 12/ 



Fdt: 



1 dH 



Fa, H^—= 



F 



F ' 



Q 



1 



dx 



2m \dt 



1 dQ 
HQ di 



(43) 



(44) 



(45) 



(46) 



(47) 



-Sf/2 



Q 
F 



? _ Sq - Sf 
" 1 + Sf/2 



l + dF/2 2H{1 + Sf/2)^' 

Since e ~ e + (5_f/2 and 6g ~ Sq — 6f at linear order in slow-roll, we find that Eqs. ([39)) . pO|) . and (|4T|) reduce to 



(48) 



ris - 1 
r 



~2i, 



16 



16e. 



(49) 
(50) 

(51) 



In the last equality of Eq. ([51]) we have used the relation e = Q/M^^, which follows from the background equation 

dH/di = -{dx/diy/{2M^i). The results ([49])-([5l]) coincide with those derived in the Einstein frame jUll^l. This 
equivalence is a consequence of the fact that both the scalar and tensor spectra are unchanged under the conformal 
transformation {Vg = 'Ps and Vt = 'Pt) [Hj. 
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Under the slow-roll conditions {\<Pxldt^\ ^ \iHdx/ dt\ and [dxl dt)"^ jl ^ [/) the background equations are ap- 



proximately given by 



We then have 



dt 



X ■ 



Q M^i fu, 



^4 



u 



pi 



Ex 
U 



u 



XX 



U 



The observables (P ^ - ([5T|) can be explicitly written as 

\,,,^,lLxx^^ 



ris - 1 



-3 A/, 



pi 



U 



pi" 



U 



V, 



V 



1/2 F 



V F 



' B 



V 



U 



(52) 
(53) 

(54) 
(55) 



In the Jordan frame the number of e-foldings from the time t (with the field value (/>) to the time at the end of 
inflation (with the field value 0/) is given by 



N 



Hdt 



- - 1 F 
Hdt + - In — 
2 Ff 



(56) 



where Ff = F{(j)f). Note that in the last equality we have used Eq. The scales relevant to the CMB temperature 
anisotropics correspond to = 50-60 [s^ . The number of e-foldings in the Einstein frame should be equivalent to 
that in the Jordan frame by properly choosing some reference length scale [40| . Using the slow- roll approximation in 
the Einstein frame, the frame-independent quantity ((56)) can be written as 



N : 



U 



xf K^u^ 



dx 



2 Ff 



(57) 



which we will use in the following sections. 



V. INFLATION WITH NONMINIMAL COUPLING AND FIELD COUPLING WITH THE KINETIC 

TERM 



In this section we study, in turn, models with the nonminimal coupling C,(j)'^R/2 and the non-canonical kinetic term 
ijj{4))X. These models are described by the action 



M. 



1 



In this case the function F is given by 



F = 1-Cx^ 



<f>/Mpi . 



(58) 



(59) 



Note that in our notation the conformal coupling corresponds to C = 1/6. For the canonical field with w(0) — 1, 
the observational constraints were studied for the chaotic potential of the type (jS]) by using the WMAP lyr data 
combined with the large-scale structure data [23| • Recently the observational compatibility of this type of potential as 
well as V{(j)) ~ X{(jP — v'^)'^ /4 was examined in Ref. [2^ by using the WMAP 7yr data. The latter potential appears 
in the context of Higgs inflation with the electroweak scale v ^ 10'^ GeV If the nonminimal coupling is negative 
with Id ^ 1, it is possible to use the Higgs field as an inflaton because the self coupling A can be of the order of 
0.01-0.1 from the WMAP normalization. Since the field 4> is much larger than the electroweak scale during inflation, 
the observational prediction of the potential V{(l)) = \{4>^ — v^Y /A is very similar to that of the potential ([3]) with 
p = 4. 

In this work we shall take into account the non-canonical kinetic term Lij{(j))X in addition to the nonminimal coupling 
C,4>^R/2. We provide general formulae for rig, r. and rit in terms of the function of a; = 4>/M-p\ and then apply them 
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to the cases where = constant and where the exponential couphng = e'^'^/^p' is present. In the Einstein 
frame this potential takes the form 

^ = (1 - Qx^y ■ 



For p < A this has a local maximum at x = ■\/p/[(4 — p)\C\\ and hence the nonminimal coupling makes it more difficult 
to realise inflation. If p = 4 the potential ([50)1 is asymptotically flat in the region (f> ^ Mp\. If p > 4 the potential 
does not possess a local maximum, but for p > 5 + vTS inflation docs not occur. 
From Eqs. and ([55]) it follows that 

~ .^^ ^ J (p-4)'(6C-cj)(Ca:')^ + (24cj-14pcj + 3p=^cj + 24K-12/C)(Ca:')^ 



uj + (6C — uj)(lx'^Yx'^ 

+ (-8w + 4pw- 3p2a; + 24pC + 6p\)Ca;^ + + 2) - ^wa;(l - Ca;^)^[(p - 4)Cx^ "P] , (61) 



where 



8b+(4-p)Cx^]^ 

'^"*-x2[^+(6C-u;)Ca;2]' ^^^^ 



= — . (63) 



For the dilatonic coupling (jj{4>) = e^'^^^'-^^ the parameter ^ is constant. Using the approximate equations (|52p . the 
scalar power spectrum is given by 

^ ^ xP+^[Qex^+^{^^Cx^)] 

' ~ 127r2Afp6 127r2M4 (1 - (a^^)^ [p + (4 - p^x'^Y ' 

The WMAP normalization corresponds to Vs — 2.4 x 10~^ at the scale k = 0.002 Mpc~^. In the following we shall 
flrst consider the nonminimally coupled theories with = and then proceed to the case in which the dilatonic kinetic 
term e^'^/*^p'X is present. 



A. Effect of the nonminimal coupling (^<^^i?/2 with constant cj 



We flrst discuss the effect of the nonminimal coupling for the theories with 



(65) 



in which case w is constant. Introducing a new fleld (p = y/ujcf) the kinetic term loX reduces to the canonical form 
— g'"^9^iy99iyiy9/2. Thcu the nonminimal coupling C,4P'R/2 can be written as C,'f'^R/2, where C = C/"^- The potential 
V^tp) = Vb(0/Mpi)^' for the scalar fleld ip, takes the power-law form V = Vo{(p / M^xY , where Vb = Vo/w^^^. This means 
that these theories reduce to nonminimally coupled theories with a; = 1 in terms of the fleld Lp. The ratio C, = C,/uj 
characterizes the effect of the nonminimal coupling on the inflationary observables ris, ?T.t, and r, and Vb = Vb/w^^'^ 
sets the scale for the scalar power spectrum. 

From Eq. (|57p the number of c-foldings is given by 



TV 



N 



J-6C 



ip- 



{x^ 



A)Cx) 




3p< 

P- 


-2oj 

-4 




1 - 


Cx^ 


4)C.t2 


-p 






1 - 


C4- 


4)- 


4 


1 






{p-- 


= 4), 


1 


-c4 





(P ^ 4) , 



where Xf = 4>f/M-p\. The result (|67l) can be also reproduced by taking the limit p 
end of inflation by the condition e = 1, which gives 



4 in Eq. 



2 _ w - Cp(4 - p) - ^{uo-2pO{lo-&pO 



C[C(4-p)2 + 2(c.-6C)] 



(66) 
(67) 

We identify the 
(68) 
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Let us consider the limits where ^ 1- We imphcitly assume that w is not different from the order of 1. We 

expand the right hand side of Eqs. (j66p and (j67p up to first order in C, and then solve them for x by using Eq. 
This gives 



p{p + 4iV) 
2u 



1 



8(p-4)jV^ + 4p(p~6)A^+p^(p-8) C 
2{p + AN) w 



(69) 



which is valid for both p ^ A and p = 4. The spectral index ((6T|) and the tensor-to-scalar ratio ((62)) are approximately 
given by 



n.s - 1 



2(p + 2) 
p + 4iV 

16p 



1 



4(p-2)(p- 12)Ar2 + 2p(p2 - 12p -f 28)Ar + ^^(12 - p) C 



p + AN 



{p + AN){p + 2) 

2N{2{p - 12)N + p{p - 10)) C ' 
p-\-AN uj 



(70) 
(71) 



This shows that the effect of the nonminimal coupling appears in terms of the ratio C,/uj. 
Substituting Eq. (|69| into Eq. (|64|) and expanding it up to first order in C, it follows that 



Vq p + AN 
M^j 247r2p 

2.4 X 10"^, 



p{p + 4iV) 



p/2 



1 - 



p4 + [AN - 12)p3 + (87V2 - 64iV)p2 + (80iV - 112N^)p + 192iV2 

4(p + 4Ar) u; 



(72) 



around N = 55. In the absence of the nonminimal coupling (C = 0) the WMAP normalization for the canonical scalar 
field if gives m ~ 6.8 x IQ-'^Mpi for p = 2 (where % = m'^M'^j2) and A ~ 2.0 x lO'^^ for p = 4 (where Vb = XM^JA). 
If ^ 0, then the inside of the last parenthesis in Eq. ([7^ is approximately given by l+AC,/oj for p = 2 and 1 -|-460C/w 
for p ~ A. As long as \C,/lij\ ^ 1, the order of Vg is not subject to change by the presence of the nonminimal coupling. 

In the following we derive the numerical values of rig and r for p = 2 and p = A separately to compare the models 
with observations. 



In order to obtain the theoretical values of ris and r for p = 2, we numerically solve the background equations of 
motion in the Jordan frame by identifying the end of inflation under the condition (pS]) . We derive the numerical 
values of x corresponding to the number of e-foldings iV = 55 and then evaluate Ug and r by using the formulas (|6ip 
and 

In Fig. [1] we show the la and 2a observational contours constrained by the joint data analysis of WMAP 7yr 
0, Baryon Acoustic Oscillations (BAO) |4l|, and the Hubble constant measurement (HST) [l^. This is derived by 
varying the two parameters and r with the consistency relation r = — 8nt [see Eq. (|4ip ]. Since the runnings of 
scalar and tensor spectral indices are suppressed to be of the order of e^, they are set to be in the likelihood analysis. 
These results are valid for the theories with ^ = = G. 

In the limit |C| < 1, Eqs. ^ and dTT]) give 



2A^ + 

16 

2iV + 1 



1 - 



4iV + 5 C 2 (l04iV4 ^ i60iV3 + MN"^ - 30 TV - 9) 



1 



2N + lu 
AN{bN + A) C 
2nTi w 



3(2iV 



(73) 
(74) 



For the scalar index we have included the second-order correction in C,/ijj because the dominant contribution to the 
first-order term in C/w in Eq. (|70p vanishes for p = 2. In the absence of the nonminimal coupling (^ = 0) one has 
Us = 0.964 and r = 0.144 for TV = 55, which is inside the 2a observational bound (sec Fig. [1]). A positive nonminimal 
coupling leads to an increase of r relative to the case C = 0. Since r is bounded from above observationally, this puts 
an upper bound on the positive value of C. The negative nonminimal coupling gives rise to the deviation from the 
scale-invariant spectrum (rtg = 1) and the decrease of r. 

From the observational constraints on rig we can place the bound on the negative nonminimal coupling. We find 
that the ratio C/w is constrained to be 



7.0 X 10"'"' < C/w < 7.0 X 10' 



(95 % CL) , 



(75) 
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Figure 1: la and 2a observational contours in the (ns,r) plane constrained by the joint data analysis of WMAP 7yr, BAO, 
and HST with the pivot scale fco = 0.002 Mpc~^ . Shown also are the theoretical predictions for the potential V{(j)) — m'^(f>'^/2 
for A'^ = 55 in three cases: (a) constant to (i.e. /i = 0) in the presence of the nonminimal coupling with ^/uj = 

0.001, 0, -0.001, -0.005, -0.01 (from top to bottom), (b) the exponential coupling e'"^/^^p'X with fi = -0.05, 0, 0.1, 1, 10 (from 
top to bottom) in the absence of the nonminimal coupling, and (c) the exponential coupling e'^'^^^p'Jt (i.e. = 1) in the 
presence of the nonminimal coupling with ^ — 0.03, 0.01, —0.05, —0.1 (from top to bottom). 



which agrees with that derived in Ref. [23 for = 1. The lower bound in Eq. ([75]) is slightly tighter than the 
constraint ( > —1.1 x 10^^ (with cj = 1) [27| obtained by using the WMAP lyr data combined with the large-scale 
structure data. 



We proceed to the case of the self-coupling inflaton potential V{4)) ~ A0*/4. In the regime IC/'^I ^ 1 Eqs. (iTOl) and 
([7T|) give 



- 1 



iV+ 1 
16 



N + 1 



1 



4{2N^ + N-4:) 
^ 3(iV+ 1) w 

4iV(27V + 3) C 
^ N + 1 w 



(76) 
(77) 



In the absence of the nonminimal coupling one has rig = 0.946 and r ~ 0.286 for N = 55, which is outside the 
2cr observational bound (see Fig. [2]). The presence of the negative nonminimal coupling leads to the increase of rig, 
whereas r gets smaller. Hence it is possible for the self-coupling inflaton potential to be consistent with observations. 
From the joint data analysis of WMAP 7yr, BAO, and HST the nonminimal coupling is constrained to be 



C/uj < -2.0 X 10 



-3 



(95 % CL) , 



(78) 



which is tighter than the bound C < —3.0 x 10^* (with ui = 1) derived in Ref. |27| . 

In another limit where IC/"^! ^ oo, inflation is realised by the flat potential U in the Einstein frame in the regime 
X > 1. In this case one has xj ~ -2^/3/(3C) and N ~ -3Ca;^/4 from Eqs. dBH]) and ^T^i, respectively. From Eqs. (HI]) 
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Figure 2: The same observational constraints as shown in Fig.[Tl with the theoretical prediction of the potential V{(l>) — \(f)* /'i 
for N = 55. Each curve corresponds to (a) constant uj (i.e. ^ = 0) in the presence of the nonminimal coupling C,(j>'^R/2 with 
C/w = 0, -0.001, -0.005, -0.01, -0.03 (from top to bottom), (b) the exponential coupling e'"*''^"'i-' X with = 0.1, 1, 10 (from 
top to bottom) in the absence of the nonminimal coupling, and (c) the exponential coupling e'^^*^p'X (i.e. /i = 1) in the presence 
of the nonminimal coupling with = 0.03, 0.02, —0.03 (from top to bottom). The label "FU" corresponds to the Fakir-Unruh 
scenario [31 with ( ^ —oo. 



and (|62p the leading contributions to Ug and r in the regime iV 3> 1 are 

n,-l ~ -2/N, (79) 
r ~ 12/N'^. (80) 

As long as \(\ is sufficiently large relative to w, the effect of the term cj appears only as the next order corrections to 
([75)1 and dBHI) with the order of uj/{CN^). For N = 55 one has ?is = 0.964 and r = 0.004 from ^ and dSD]), which 
are well inside the la observational bound. 

In the regime \C/Ld\ ^ 1 the power spectrum ([M)) reduces to Vs — AAf^/(727r^(^^), so that the WMAP normalization 
Vs ^ 2.4 X 10"^ at TV = 55 gives 

A/C' ^ 5.6 X 10-1°. (81) 

For large negative nonminimal couplings, such as C ^ —10^, the self coupling A can be of the order of 10^^. This 
property was used in the context of Higgs inflation. 



B. Effect of the non-canonical kinetic term e^^^^'^'^'X with C = 

Let us consider the case in which the field-dependent coupling uj{(t)) = e^"^/^^?' with the kinetic energy X is present, 
without taking into account the nonminimal coupling (C = 0). After the field settles down to the potential minimum, 
(/) = 0, the coupling w(0) — )• 1 and one recovers the standard kinetic energy X. 

The number of e-foldings (j57p is given by 

^ = J- [(^a; - l)e^^ - Uxf - l)e^^-f] . (82) 
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Since e = / (2x^0;), we can estimate s/ by setting e = 1: 



2VF(y2|/i|p/4) 



(83) 



where W is the Lambert's function [4^- The scalar spectral index (|5T|) and the tensor-to-scalar ratio ([S^ are 



ns - 1 = 



0^2 



(p + 2 + ^x) , 



In the limit |^| <C 1, using Eq. (|83|. one can rewrite Eq. (|82|) in the form 



N ■ 



2a;2 - p2 



4p 24p 
which can be solved for X, cLS 

2 

By replacing this relation into Eqs. (|M)) and ([55)) . we obtain 

p + 2 



(84) 
(85) 

(86) 

(87) 



ris - 1 



2iV 



r ~ — 1 

N 



3(p-f2) 

3 ; ' 



(89) 



which are valid up to the first order in fi. The presence of the positive /x leads to the approach to the scale-invariant 
spectrum, whereas r gets smaller. In Figs. [T] and [2] we plot the theoretical values of ris and r in the {ns,r) plane for 
p = 2 and p = 4, respectively, with several different values of ^. These are derived numerically by integrating the 
background equations without using the approximation given above (because the approximation loses its validity for 
/X > 1). Interestingly the models with large positive values of fi can be favoured observationally. On the other hand, 
the models with negative n lead to the deviation from the observationally allowed region. The joint observational 
constraints from WMAP 7yr, BAO, and HST give the following bounds on fi: 



fi > -0.04 (95% CL) for p 2 , 
Ai>0.2 (95% CL) for p^A. 

Let us now consider the limit where fi^ 1. In this regime the condition fix ^ 1 is satisfied, so that N ~ xe^ 
from Eq. (|82p . Then the scalar index and the tensor-to-scalar ratio (|55|) reduce to 



1 



1 

8p 1 

'nJIx 



(90) 
(91) 



(92) 
(93) 



For a given A'', /ix increases for larger /i. This means that, in the limit ^ ^ 1, ng and r approach Us ^ 1 — 1/N ~ 0.982 
(for A'' = 55) and r — >■ 0, respectively, which is inside the la observational bound. We have also confirmed numerically 
that inflation is followed by a reheating phase with oscillations of (j). 



C. Combined effects of the nonminimal coupling C,(t)'^R/2 and the non-canonical kinetic term e^'^^'^'i'^X 

Finally we consider the case in which both the nonminimal coupling C,(j)'^R/2 and the non-canonical kinetic term 
gM^/J^^pi taken into account. Since it is difficult to derive an analytic form for the number of e-foldings A^, we 

solve the background equations numerically to identify the values x corresponding to A^ = 55 before x = Xf which 
happens when e = 1. We then use the formulae ([54]) and (fSS]) to evaluate Us and r for given values of p, fj., and (. 
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In Fig. [T] we plot the numerical values of Us and r in the two-dimensional plane for p = 2 and fi = 1 with ( = 
0.03, 0.01, —0.05, —0.1. The presence of the term e''"^/*^p'X with /i > leads to the compatibility of the nonminimally 
coupled models with larger values of |C| than those for /i = and = 1. When /i = 1 we find that the nonminimal 
coupling is constrained to be 

- 0.12 < C < 0.035 (95% CL) , (94) 

which is wider than the range ()75|) . For values of |C| larger than the bounds given by ()94|) the effect of the nonminimal 
coupling is more important than that of the non-canonical kinetic term. 

For p = 4 and fi ~ a positive nonminimal coupling is not allowed observationally because both {ng— 1| and r tend 
to be larger than those for C = 0. However, the non-canonical kinetic term with /i > allows the compatibility of the 
positive nonminimally coupled models with observations (sec Fig. [5]). If /i = 1, C is constrained to be 

C < 0.025 (95% CL) . (95) 

For fJ- ^ 1 the models with C < are within the la observational bound. In the limit of the largely negative nonminimal 
coupling (Id ^ 1), the scalar index and the tensor-to-scalar ratio are given by Eqs. ([79| and ((80|). 

For fi larger than the order of 1 the effect of the non-canonical kinetic term tends to be more important. In the 
limit that fi^ 1 with a finite value of C, (where |C| S l)i "-s and r approach the values given by Eqs. and (^5)) . 



VI. INFLATION IN THE CONTEXT OF BRANS-DICKE THEORIES 



Let us proceed to Brans-Dickc (BD) theory [4J| with the action 



S = / d^x^ 



].M^,c^R+^Lo^T,X-V{<jy) 



(96) 



where wbd is the BD parameter. Here we have introduced the reduced Planck mass M-p\ in the first two terms, so 
that the field has a dimension of mass. Under the conformal transformation we obtain the action in the 
Einstein frame with 

F=-^ = e'^'^/^^p' , U = e-2^>^/*^pi V , (97) 
Mpi 

where 

1/V3/2 + CJBD. (98) 
The integration constant for the field x is chosen such that x = corresponds to cj) = Mpi. 



A. Case of the power-law potential 

Let us consider the power-law potential ([3]) in the Jordan frame. In the Einstein frame the potential is given by 

U{x) = Voe^^/^'"^ , A ^ J- '^ . (99) 

y/3/2 + ujbd 

For the BD parameter of the order of 1 infiation does not occur unless p is close to 2. However, for wbd ^ 1, it is 
possible to realise |A| ^ 1 even Hp is away from 2. 

On using Eqs. ([54|) and ((55|) for the potential ((99|). it follows that 



Hs-l = -X\ (100) 

r = -8nt = SA^. (101) 

The CMB likelihood analysis using the data of WMAP 7yr Q combined with BAO [U and HST gives the 
following bound on A [i^ : 

0.09 < A < 0.23 (95% CL). (102) 



14 



This translates into the constraint on ujbd' 

19{p - 2f - 3/2 < CJBD < 123(p - 2)2 - 3/2 . (103) 

The reason why the p = 2 case [Le. A — 0) is disfavoured is that the Harrison- Zel'dovich spectrum {ug = 1 and r = 0) 
is in tension with observations Q- If p = 4, Eq. p03p gives the bound 75 < ojbd < 491. 

The exponential potential in the Einstein frame does not lead to the end of inflation, so the above scenario has 
to be modified in a way that the potential has a minimum to lead to a successful reheating. In the following we 
shall consider the modification of the power-law potential in the Jordan frame, such that inflation ends as in the 
Starobinsky's /(i?) model [l|. 



B. Models including the Starobinsky's /(i?) scenario 
The f{R) theory with the action 

/j\,/2 

is equivalent to BD theory with ojbd = [i^. In fact the action (|104p can be written as 



S 



I- 



Ml, 



F{<j>)R ~ F(0) 



where 



df 



M, 



pi 



dR 



Af 



dR 



In Starobinsky's model f{R) = R + R'^/{6A'P), we have R = 3AP{(p/Mpi - 1) and 



F(0) 



3AP 



-A/pi)^ 



We consider the following more general potential 

F(<^) = Fo(0-Afpi)^ 



(104) 



(105) 



(106) 



(107) 



(108) 



with arbitrary values of wrd, so that the model f{R) = R + R^/{6M^) is covered as a special case with p = 2 and 
■^BD = 0. The potential in the Einstein frame reads 



(109) 



where /i is defined in Eq. (|98p . For [wbdI ~ i.e. jj. ^ 0(1), inflation occurs in the regime x ^ -^^pi (including 

the Starobinsky's f{R) model). The behaviour of the potential (|109p depends on the values of p: 

• If p = 2 the potential (|109p becomes constant for x ^ -^'^pi- We note that it is not necessary for p to exactly 
equal 2 for this behaviour to occur. Since U is approximated a,s U oc x^ in. the regime x Mpi, inflation is 
followed by a successful reheating. 

• For p > 2 the held rolls down the potential towards x = 0. 

• When p < 2 the held rolls down the potential towards x ~ +oo or towards x = 0. In the latter case the potential 
does not have a minimum at = 0. As a result, reheating is problematic for p < 2. 



From Eqs. and it follows that 



Hs- 1 = - 



[4 + 2{3p - 4)F +{p- 2)2^2] 
8n^[2 + ip~2)F]^ 



(110) 
(111) 
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The number of e-foldings (|57p reads 



»= ^(f-F,) + i(l--i)l„(^) (p = 2). (112) 



where Ff is the value of F at the end of inflation. Using the criterion e = 1 for the end of inflation, we have 

Ff^ ^ + ^. (114) 



i. Case: p — 2 

Let us consider the case p = 2. For the theories with |wbd| ^ C'(l) (i-C- M ^ C'(l)) one has = 1 + \/2a' = C'(l) 
and TV ~ F/{2fi'^), which means that 1 for N ^ 1. From Eqs. (jllOp and (jllll) it follows that 



F ^ TV ' 
32^2 ^ 8 _ 4(3 + 2a;BD) 
~F2~ ~ J!^ ~ TV2 



(115) 
(116) 



which are valid for —3/2 < wbd < £>(!)■ The metric f{R) gravity corresponds to ujbd = 0, which gives r ~ 12/iV^. 
This result matches with the one derived in other papers [431 ■ In the limit that wbd —3/2 the tensor-to-scalar 
ratio vanishes. The BD parameter wbd = —3/2 corresponds to Palatini f{R) gravity [HI [HI, in which case a separate 
analysis is required as in Ref. jisl ]. 

If wbd 3> 1, then one has /i ^ 1 and hence F is close to 1 even during inflation. The end of inflation is characterized 
by the condition e = 1, which gives Ff ^ 1 + %/2/i ~ 1. Then the number of e-foldings (jll2p is approximately given 
by ~ (x/Mpi)^/4. The scalar spectral index and the tensor-to-scalar ratio are 



r.32-fi.-, (118) 

which match with those for the chaotic inflation model with the potential U{(j)) = m^0^/2 [49| . 

The tensor-to-scalar ratio depends on the BD parameter wbDi while the scalar index is practically independent of 
LdBD- In Fig- El we plot the theoretical predictions of Us and r for several different values of ljbd by fixing N = 55. 
Shown also are the la and 2a observational contours constrained by the joint data analysis of WMAP 7yr Q, BAO 
HH, and HST [42|. The f{R) model f{R) = R + R^/{6M^), which corresponds to wbd = 0, is well within the la 
observational contour. While the present observations allow the large BD parameter with cjbd ^ 1, it will be of 
interest to see how the PLANCK satellite [lOl can provide an upper bound on wbd- 

Using the approximate relation F ~ 2ii^N + Ff following from Eq. (|112p . the WMAP normalization for the scalar 
power spectrum Vs = / {12n^ M^^U"^^) is given by 

-24x10'^ (119) 

around 55. Since wbd = and l/\/3/2 for the f{R) model /(i?) = R + R'^/{6M^), the mass scale M is 
constrained to be Af ~ 3 x 10^^ GeV. The energy scale Vq is different depending on the BD parameter. 



2. Case: p / 2 



We proceed to the case p ^ 2. For the BD parameter wbd of the order of unity, the number of e-foldings (|113p 
cannot be much greater than 1 unless F is enormously larger than Ff (^ 0{1)). If ^ 1, then Eqs. (jllOp and (|llip 
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Figure 3: Icr and 2a observational contours in the (ris, r) plane constrained by the joint data analysis of WMAP 7yr, BAO, and 
HST with the pivot scale fco = 0.002 Mpc~^ (logarithmic scale for the vertical line). The dotted points show the theoretical 
predictions for the BD theories with the potential V{(j]) = Vb(<^ — Mpi)'^. The number of e-foldings is chosen to be A'^ = 55. 
From bottom to top the points correspond to ij;bd = — 1.4, -1, 0, 10, 10^ 10^ 10" and ujbo — oo, where wbd = represents the 
model f(R) = R + B? / {GM'^). For larger i^bd the two observables Us and r approach those for the chaotic inflation with the 
quadratic potential m?(jP' /2. 



give 

ris-l ~ V(p-2)^ (120) 
r ~ 8/i2(p-2)2. (121) 

Since /i ^ 0{1) the results (|120p and (|12ip mean that for smaU wbd both the scalar index and the tcnsor-to-scalar 
ratio are incompatible with observations apart from the case where p is close to 2. On reflection this is unsurprising 
since only for p w 2 is there a flat region of the potential that will give slow-roll along with its signatures of near scale 
invariance and suppressed tensor modes. 

When wbd S> 1 one has /i ^ 1 and hence F = e^'^/^^p' is close to 1. Then Eqs. (jllOp and (jllip give 

-p(p+2)-fi^-^, (122) 



X2 2N 

(123) 



where we have used the approximate relation N ~ / i'^P^'Ipi)- These results match with those of chaotic inflation 
with the potential (jSj. The case p = 4 is excluded observationally both in the regimes wbd ^ 1 and ojbd = C'(l)- 
Even for other values of ojbd it is difflcult to satisfy observational constraints unless p is close to 2. 
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VII. INFLATION IN THE PRESENCE OF A GAUSS-BONNET TERM 



In this section we study the effects of the Gauss-Bonnet (GB) term on the chaotic inflationary scenario, described 
by the action 



S = J d'^x^/^ 



Ml 



R + x -v((b)-^{(f,)g 



(124) 



In order to confront the model with observations, it is convenient to rewrite inflationary observables in terms of the 
following slow-roll parameters: 



V 



V 



The background equations are 

3M^jjf2 ^ 02^2 + y + 2m^i , 

At linear order Eqs. (pSj) and give 

Cs = (5x , e = + 4(5^ . 
From Eqs. (|126p and (|127p the potential V and its derivative can be expressed as 



V = SM^H^ ( 1 - ^e., - 



1 s 

3-e+-??s + 12^(l-e) 

Z Cc 



= -Hc^ 

Taking the leading-order contribution in Eq. (|130p . it follows that 



which lead to 



e., 1 



-Tjs - 2e^ - 16(55 ^ 8(5^ + -rj^ - ?/e 



;Vs 1 



4(5, 



2e, + 



4+ -(8(5^-775) 



From this we obtain the inversion formulas 

1 

e, ~ — 
2 



, - 8(5^ + ^e? - i6e,,(54 , 

2 {t?^ - 2e, - 4(5^ [4 + (8i5g - 7?g)/e,]} 
1 - 4.S^/es 



where we have taken the positive sign in Eq. (|135p to reproduce — > e^, for (5^ — > 0. 
From Eqs. (p8|) . (|3ip . and ((32)) the inflationary observables are given by 

ris - 1 = -2es - Vs - 8(5^ , 
rit = -2es - 8^5 , 
r = 16es , 



(125) 

(126) 
(127) 

(128) 

(129) 
(130) 

(131) 
(132) 

(133) 
(134) 

(135) 
(136) 



(137) 
(138) 
(139) 
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which are written in terms of the four variables: ev^ ??v, 1^5, and 77^. By specifying the functional forms of and 
^(0), we can reduce the number of those variables. For the chaotic inflation potential ^ one has tv = (p^/2)(Mpi/0)^ 
and r\v = vip ^ Ijl-^pi/*/")^; so that they are related with each other via the relation 



fcll... (140) 



For the Gauss-Bonnet coupling, we take 



m = i^e^^l^'-' , (141) 
where ^0 and fi are constants. It then follows that 

77. = -2e, + - 8S^ ± ^V2il , (142) 

where the plus and minus signs correspond to > and < 0, respectively. Combining Eq. (|142p with Eq. (|136p . 
we obtain 



Vs - -2 



77v - 2es + 4:6^ ^ - 6 



(143) 



Substituting Eq. (|143p into Eq. (|137p and choosing the negative sign of 0, the scalar spectral index can be written as 

7^, - 1 ^ _6e, + 27y, - 8^, ^7 + [^^^^±1^^ , (I44) 

where 

2(p-l) 2(p-l) / 45. V , 
P P \ <^s J 

For fixed values of p and fi one can carry out the CMB likelihood analysis in terms of Ug, r, and nt by varying the 
two parameters Cs and S^. 

In Fig. |4]the observational constraints on the parameters and = S^/eg are plotted for p ~ 2 and /i = 1. We 
run the Cosmological Monte Carlo (CosmoMC) code [5l| with the data of WMAP 7yr Q combined with large-scale 
structure Q (including BAG [U), HST [42|, Supernovae type la (SN la) [HI, and Big Bang Nucleosynthesis (BBN) 
[ssl j. by assuming a AGDM universe. The ratio is constrained to be \r^\ < 0.1 (95% GL), which means that the 
effect of the GB term needs to be suppressed. Hence the energy scale Vq is similar to that in the standard chaotic 
inflation. 

From Fig. |4]we find that the slow-roll parameter eg is bounded to be < 0.025 (95 % GL). In the presence of the 
GB term the small values of can give rise to the scalar index close to 77,^ = 0.96. For example, when = 0.002, 
~ 0.05, /i = 1, and p = 2, one has 77,5 = 0.962 from Eq. (|144p . This is different from the standard chaotic 
inflation in which the small values of lead to the spectrum close to the Harrison- Zel'dovich one (which is not 
favored observationally) . Hence the allowed range of Eg tends to be wider in the presence of the GB coupling. 

Let us estimate the two observables 7t,s and r in terms of the number of e- foldings N under the condition \r^ \ 1. 
Since ~ - 8S^ from Eq. Eqs. and reduce to 

ns~l ^ - ( 2 + - ) - 8S^ (1 + IJ\ —] , (146) 



PJ V V e 

r ~ 16e,. (^1 - . (147) 

From Eqs. and p^ one has if/0 ~ -(1 + 4:6^/es)V/{M^iV^). Using the relation e., = Sx and the definition 

of S^, it follows that H/cp ~ —(j)/{pM^i + 8i7^^^00). Since we are considering the case where H^^^^(j)/M^i <^ 1, the 
number of c-foldings for the potential ^ is 



N= dcj,^^-^+Np, where Np = -^^ e^^x^+Hx. (148) 

' ' 2p 3p^ M^j 
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Figure 4: la (inside) and 2cr (outside) observational contours in the (es,r^) plane for the potential V{(t>) — m?(j)'^ /2 with /i = 1 
(r^ = 5s^/e,). We use the data of WMAP 7 yr, LSS (including BAO), HST, SN la, and BBN with the pivot scale fco = 0.002 
Mpc-^ 



Here a;/ is the value of a; = 0/A/pi at the end of inflation. We identify the end of inflation by the condition ev = 1, 
i.e. Xf ^ pIV^- When p ~ 2 and p = 4, Eq. (|148p is integrated to give 



N2 



{e''^[/^4a;4 4- ^{-[i^x"^ S/i^^ - %[ix + 6)] - 4e^^ [m^* - 2^2/.*^ + 6(m^ - + 1)]} , (149) 



iV4 = <! e^'^(//'^a;'^ - 6//^x^ + 30/a;^ - 120^^a:^ + Zm^^x^ - 72Q^ix + 720) 



-16e 



32^i^ - 48%/2^^ + UOfi^in -V2) + 90/1(2^ - \/2) + 45 



(150) 



where we have set Vo = '7i^Mp[/2 for p = 2 and Vq = AA/pj/4 for p = 4. 

For positive ^ of the order of unity, the dominant contributions to Np come from the first terms in Eqs. (|149p and 
(fTSO)) . i.e. N2 ~ -m2^oa;''e^^/(3M2i) and iV4 ~ -ACoa;^e^^/24, for the scales relevant to CMB {^ix > 1). In this case 
the number of e-foldings (|148p is approximately given by 



N 



N 



!--( — ) Cx' 
3 \Mpi' 



1 



(151) 
(152) 



Since ~ — fiVop / {3M*^) x^ and = P^/(2a;^), one can express the scalar index (|146p and the tensor-to-scalar 
ratio (|147p in terms of x. By treating the ^-dependent terms in Eqs. (|15ip and (|152p as small corrections, ng and r 
can be written in terms of N: 



n„ - 1 



2 






1 


TV 




3 






1 


iV 





f' 

256 
~9~ 



m 



pi 



N 

16 

TV 



(p = 2) , (153) 
(p = 4) , (154) 
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which are vahd for positive ^ of the order of unity. If ^ > (i.e. i5j < for (j) < 0), then the effect of the GB coupHng 
leads to the approach to the scale-invariant spectrum, while r gets larger. The negative values of ^ lead to the decease 
of r, but rig deviates from 1. Since eg is approximately given by eg ~ w p'^/{8N), the scales relevant to the CMB 
anisotropics (TV = 50-60) correspond to 0.008 < < 0.01 for p = 2. Figure 2] shows that the ratio is constrained 
to be —0.04 < < 0.03 (95 % CL) for this range of e^. The self-coupling potential V{<f>) = A(/)'*/4 is not saved by 
taking into account the GB term with positive fi, because the GB coupling does not lead to the increase of ng and 
the decrease of r simultaneously. 

For negative p with — 0(1) the exponential term e^^' in Eqs. (|149p and (|150p is much smaller than 1 for the 
scales relevant to CMB {x ^ 1). In this case wc have 



N4 



4m^ ^oe 



2^^3 + 6(/i2 - \/2^ + 1) 



32/i^ - 48^2//^ + UOp^iiJ -V2)+ 90/i(2/i - V2) 4- 45 



(155) 
(156) 



The scalar index and the tensor-to-scalar ratio are approximately given by 



Us - I 

n.= - 1 



1 (i+^^i^m 

N \ N 



3 / N4-I 



N \ N 



16 , ^4-1 



(P = 4), 



(157) 
(158) 



where we have assumed N ^ Np, and ignored the exponential term e^^. When fj, < 0, one can show that N2 and N4 
in Eqs. (|155p and (|156p are positive for ^0 > and negative for ^0 < 0. In the latter case the presence of the GB term 
leads to the approach to the Harrison Zel'dovich spectrum. In fact, such a scenario was discussed in Ref. [2l|. Since 
m/Mgi and A are much smaller than 1 by the WMAP normalization {m/Mpi ~ 6.8 x 10~^ and A ~ 2.0 x 10"-^"^), one 
has \Np\ ^ 1 for |^o| smaller than the order of 1. For /x < the effect of the GB term on the inflationary observables 
appears for very large values of ^0 such as |^o| ^ 10^° [2l| . 

VIII. G-INFLATION WITH A FIELD POTENTIAL 

Finally we study chaotic inflation in the presence of the Galileon-like self-interaction G((f>, X)D(f> (called "G- 
inflation"). We specify the functional form of G{(j>,X), as 



G(0,X)==$(0)X", $(0) 



(159) 



where n and /x are constants, and 6 = ±1. The constant M has a dimension of mass with AI > 0. Here we have 
introduced the exponential form for $ motivated by the dilaton coupling in the low-energy effective bosonic string 
theory. Wc also consider the power-law function by generalizing previous studies [30j . 
Equations ([5]) and (jT]) can be written as 



V = 3Mp2 [l--5x 



3i70 n + (3 - e) 



2(5gx 

6gx _ 
Sx 



3n Ox 



1 + 6n 



'GX 



2{n+l) 



^G4> 



X 



(160) 
(161) 



To compare with observations, we seek an expression for rig — 1 = Vsl {HVs) in terms of a minimal set of independent 
slow-roll parameters. Since Vs = /{Stt^Qc"^), it is important to find the expressions for Q and Cg. From Eqs. ([15]) 
and (1161) it follows that 



Q_ 



Sx + GnScx -2{n+l)SG4, + ^^gx 



^GX 



Sx + 2(2 4- nS^)SGx + 2{n - 1)Sg4> - ^> 



GX 



■Jx 



6nS, 



GX 



2{n+l)SG^ + 35l^ 



(162) 
(163) 
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where we have used the relations Xgx = n — 1, and Xq^ = n. Hence, we can derive an exact expression for rig — 1 
in terms of the slow-roll parameters entering Eqs. (|162[) and (|163p and their first derivatives, which introduce other 
slow-roll parameters, however these are not all independent. We shall now discuss the relations which reduce the 
number of independent slow-roll parameters. 

For the choice of the function $ in Eq. (|159p we have 



\/2r. 



5g^ = ±^SgxVSx, (164) 



where the ± signs in this expression are compatible with those in the expression (j) ~ ±\/2Mp\Hy/ Sx ■ Equation (jl64p 
shows that 6041 is in general suppressed relative to Sgx- This relation also implies that 

VG4> = TFF — = ^ + mx , (165) 



where 



From the definition of 5x and 6gx we obtain 

r,x = 2{1 - Sgx)S^ + 2Sx + 65gx - 4Sg^ , (167) 
TjGX = i2n+l-6Gx)S4,±^iV2S^ + Sx-2SG4, + SSGx, (168) 

where we have used the relation 

e = Sx + 3Sgx - 26g^ - <5^<5gx ■ (169) 

It should be noted that the four relations (fT64| - (fT69| are ah exact. From Eqs. (fT65l) . (fT67l) . (fTeS]) with Eq. (fTM)) we 
find that t?g0, , and rjGx can be expressed in terms of the three slow- roll parameters 5^, Sx, and Sgx- 

We finally use a last constraint coming from the fact that we have chosen a power-law form ([3]) for the potential. 
The relation (|140[) between ey and r/y leads to 



P-IV4 



HV4, p HV" 

This equation can be used to set the last constraint on the slow-roll variables. At lowest order we have 



(170) 



{Sx +iSGx)[{2- p)Sx +^Sgx] 3V2(5gxv^ 2{n - \)SxSgx{Sx ~ ^Sgx) 2 , 3/2\ 

7x — r«"^^ — ^ T — rT^ — ^^ tj — — \2 ^ +c(e' )■ (i7i) 

p{Sx + bnSGx ) Sx + ouSgx n[Sx + bnScx ) 

Using this relation we can express 7700 1 Vx, and i]gx in terms of two slow-roll parameters Sx and Sgx- 

We are now ready to explicitly calculate the scalar index rig — 1 = — 2e — Sq — 3s, where Sq ~ Q/{HQ) and 
s = Cs/{Hcs) are evaluated by taking the time derivatives of Eqs. (|162p and (|163p . This gives 

"s-1 -— ^('^^+3'^^^-) [4(p + 2) + <5i^G^[(3p-6)n^ + (12p + 27)n + 4p + 8] 

p (Ox + 4 Sgx ) [Sx + 6n Sgx ) 

+ 5xSgx[{^'^P + 30)^2 -f (54p + 105)n + 6] + 72n<5^_y(3np + 2n+ 1)] 
^ 3y2fcxv^[(7n + 2) <5j^(5gx + + 2^nSlx\ 
((Jjf + 6ri(5Gx) [Sx+^Sgx) 

'^1^1^ [Sx^ + (9 n + 8 + 6 n3 - 24 n^) S\Sgx 

n{Sx + QnSGx) {^Sgx+Sx) 

+ (4 - 99 -f 54 n - 42 n^) S^S^x + (l32 ~ 282 - 24 - 132 n) SxS%x - 72n {r? - 3n + 6) S%x] M^ 

(172) 
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where, in order to derive this result, we have also included the terms of order 0(e^/^) not shown in Eq. (|171|) . Again 
the ± signs in the term involving n in Eq. (|172p are compatible with those in the expression for </> in terms of Sx- 
The tensor-to-scalar ratio ([32]) and the tensor index ([31]) are approximately given by 



r ^ 16 -jyj, (173 

nt ~ -2(^x + 3(5gx), (174) 



and the scalar propagation speed squared is 



-I - r ^."'^r ■ (175) 



If I <^GX I ^ , then these obscrvablcs reduce to 



1 - -fc^<5.±3V2nM^, (176) 



r ~ 16(5jf ~ -8nt . (177) 

On the other hand, in the limit where Sgx 3> ^x, one has 

3(3?ip + 271 + 1) , M /7- 

ns-1 ~ dcx ± Vf>x , (178) 

V2n 

64 /e"^ 32 |6" 8.7 

which agree with the results in Rcf. ^| derived for n — \ and /i = 0. 

To be concrete, in the following discussion we focus on the theories with n— 1, ^ 7^ 0, and 6* = — 1. Then ^qx > 
for < 0, so that the conditions for the avoidance of ghosts and Laplacian instabilities {Q > and > 0) are 
always satisfied. In this case we need to take the minus sign for the term fi in Eqs. (|172p . (|176p . and (|178p . Since 
~ —3H(j>{l + 377 (/)<&) from Eq. (|16ip . the field velocity corresponding to < is 



v/1 -4$V"rf,- 1 , , 

^ (180) 



6H^ 

r2 J f 2 



where we used $ < 0. Employing the approximate relation V ~ 3H AI^^, the two slow- roll parameters Sx and 6gx 
can be expressed in terms of 0, as 



3X — 



8V"2$2 ' 5 

The number of c-foldings is given by 



Sgx ^ ^iVT^^^W^ - I) ■ (181) 



Af2j ^1 - 4$v:0 - 1 7:,^ yiTiBv^^v^ - 1 ■ 



iV= / / —————d4> = 2B^ I ^ _ ^ dx, (182) 



where 



We determine the value of a;/ at the end of inflation using the condition e ~ (5^ + ZSgx = 1- 

In the limit B ^ (i.e. 5gx 0) we have e^6x^ p^/(2a;^) and TV ~ x^/{2p) - p/4, so that Eq. p77)) gives 

n.:.l-fc^, r:.^.-8n,. (184) 

In Eq. p84p we have not taken into account the contributions coming from the term /x, because we do not have 
analytic an expression for general p. Numerical calculations show that in the regime i? <C 1 both Ug and r become 
smaller for /i > 0. If < 0, then rig get smaller, whereas r increases. 
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Figure 5: Theoretical values of ris and r for the potential V{(f>) — m^(;/)^/2 in the presence of the Galileon- 
type coupling G = — (l/Af^)e'"^/*^piX with fi = 1 (solid line). The points correspond to the cases with B = 
0, 10"^/^ 10"^/'', 10"^ lO"'^'''', 10"^/^ 10-^/*,0.1, 10"^/^ 1, 10^/^ 10, 10='''^ 10^ with iV = 55. In the limit B -> oo one has 
Tis = 0.9675 and r — 0.1258. The dotted curve corresponds to the case where fj. — 0. We also show the la and 2a observational 
contours derived by the joint data analysis of WMAP 7 yr, BAO, and HST with the consistency relation r — —8nt. 
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Figure 6: Similar to Fig. (5] but for the potential V{(f>) = A(/>*/4 with /x = 1 (solid line). The points correspond to the cases with 

B = 0, lo-3/^ 10-"/*, 10-*, w-'^/\ lo-^/^ 10-"/*, lo-^ lo-"/*, lo-*/^ lo-^ lo-3/^ lo-^/*, w-^/^ lo-^/*, i, lo^/^, w^/\ 

10=*/^ 10^ with TV = 55. In the limit where B — > oo one has Us = 0.9614 and r = 0.1791. The dotted curve corresponds to the 
case where jj, = 0. The same observational contours as those in Fig. [5] are also plotted. 
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In the opposite limit where B ^ I it follows that 



(185) 



and 5x — px '^^+^^6 (25'*). In order to have « 55 for B ^ \ the integral inside the expression of N needs to be 
much smaller than 1, so that x ^ 1 for = 0(1). Using the approximation ^ 1, we have x^J'^^''^^ ~ p'^/^/ (25^) 



and 



N 



B^ 2 



.(p+3)/2 



From Eqs. ((T78)) and (|T79)) it follows that 
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2(p + 5) 



/i.T 



P + 3 



^ 3(p+l) 
(p + 3)A^ + p 

64^6 p 

9 (p + 3)iV+p 



2(p-l) 



3(p+l)(p + 5) 



fix 



P 



(186) 

(187) 
(188) 



For N — 55, in the limit where /i — > 0, one has rig = 0.9675, r — 0.1258 for p = 2 and rig — 0.9614, r — 0.1791 for 
p = 4. In the regime B ^ 1 the current observations can be consistent with both models. In the presence of the 
exponential coupling with positive /i the scalar spectral index gets larger for p > 1 , while the tensor-to-scalar ratio is 
smaller. 

In the intermediate regime between B <^ 1 and _B 3> 1 we evaluate rig and r as follows. For given values of p, ^, 
and B wc identify the field value x = (j)/M-p\ corresponding to = 55 by integrating Eq. (jl82l) numerically. We derive 
5x and 5gx from Eq. p81|) which allows us to obtain ng and r by using the formulas (|172p and (|173p . We have also 
solved the background equations numerically to the end of inflation and confirmed that the above method provides 
accurate estimation for Us and r. 

The theoretical values of Us and r for ^ = 1 are plotted in Figs. [S] and [5] (corresponding to p = 2 and p = 4, 
respectively) with several different values of B. If we choose larger B starting from i? = 0, rig decreases up to some 
value of i?, starts to increase, and flnally decreases towards the point given by Eq. p87p . Meanwhile r decreases up 
to some value of B with a minimum smaller than 0.1, before starting to increase towards the asymptotic value (jl88p . 
The above peculiar curved trajectories in the (rig, r) plane occur because of the presence of the exponential Galileon 
coupling with > 0. For fi — Q the theoretical curve can be approximated by a line that connects the two asymptotic 
points corresponding to i? and B oo, see Figs. [5] and [H 

li fi < and B is increasing, then r increases up to some value of i?, whereas rig decreases. The maximum values of 
r for /.i = —1 are about 0.35 and 0.68 for p = 2 and p = 4, respectively. If B is increased further, r starts to decrease 
towards the point given by Eq. (|188p (with rig starting to increase at some value of B). Compared to the case /i > 
this behaviour is not desirable to satisfy the observational bounds, especially for p = 4. In the following discussion 
we shall therefore focus on the case of the positive fi. 

From Eqs. (|173p and (|174p one has 



(1 + 4i?G) 



3/2 



nt 



{1 + GRay/Hl + 3Rg) 



(189) 



where Rg = Sgx/Sx- For < Rg < oo the ratio r/rit is constrained to be in the narrow range —8.71 < r/nt < —8. 
We carry out the CMB likelihood analysis in terms of Ug and r by using the two consistency relations r = —8nt and 
r = — 8.71nt. We find that the observational constraints on rig and r are similar in both cases. Hence the constraints 
using the standard consistency relation r = —8nt should be trustable even in the intermediate regime. Figure [5] shows 
that the quadratic inflaton potential is consistent with observations even in the presence of the exponential Galileon 
coupling with /i = 1. From Fig. [S] we find that the self coupling inflaton potential can also be saved by taking into 
account the exponential Galileon coupling. 

For the theoretical points shown in Figs.[5]and|6]wc can evaluate the values of Sx and Sgx corresponding to = 55. 
It is then possible to derive fitting functions that relate 5gx with Sx- The fitting function for p ~ A and /i = 1 is 
given in Eq. (jAip in the Appendix. This allows us to run the CosmoMC code in terms of one inflationary parameter 
5x- In Fig. [7] we show the 1-dimensional marginalized probability distribution for p = 4 and /i = 1 constrained by the 
joint data analysis of WMAP 7 yr, LSS, HST, SN la, and BBN. In the absence of the Gahleon coupling {5gx = 0) 
one has Sx = p/{AN + p) — 0.018 for A^ = 55, which is obscrvationally excluded. In the opposite limit of large 
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Figure 7: 1-dimensional marginalized probability distribution of the parameter Sx for G-inflation with the quartic potential 
V{4>) = A0''/4 and fj, — I. We use the fitting function ([Af |) that gives the relation between Sx and Sax for A'' = 55 in the 
regime 10~* < Sx < 0.018. The parameter Sx is constrained by the joint data analysis of WMAP 7 yr, LSS (including BAO), 
HST, SN la, and BBN with the pivot scale ko = 0.002 Mpc"^ 



Galileon coupling such that Sgx > Sx-, it follows that Sqx — p/{3[(p + 3)iV + p]} = 3.4 x 10~^ for N = 55. Since 
this case is marginally inside the 2a observational contour in Fig. [6l we find a suppressed probability distribution for 
smaller Sx in Fig. [T] The intermediate regime such as 10^"* S ^x ^ 10~^ is most favored observationally, because 
the corresponding theoretical points can be deep inside the 2a bound in Fig. |6l In Fig. |6] the theoretical point for 
B = 10"'^/^ gives Sx = 3.5 x 10~^, which actually corresponds to the highest probability in Fig. [T] Hence the effect 
of the exponential Galileon coupling works to save the self-coupling inflaton potential. 

The scalar spectrum Vs at the scale k = 0.002 Mpc~^ (for ri = 1) is subject to the WMAP normalization: 

where y = {1 + ApB^xP-'^e^^'^f/'^. In the limit that B > 1 we obtain m « GeV/M) GeV for p = 2 (with 

Vo = m^M^i/2) and A w (lO^^ GeV/A/)-^ for p 4 (with Vq = AAf^ /4), which agree with those obtained in Ref. [s^ 
for fi ^ 0. In order to have i? ^ 1 for p = 4, we require that A ^ (M/Mpi)^. Combining this with the WMAP 
normalization, the mass scale M is constrained to be M ^ 10~''Mpi. If we demand that the coupling A is smaller 
than 1, this gives another constraint M > 4 x lO^^A/pi. In the intermediate regime between B 1 and _B ^ 1 we 
need to solve Eq. (|190|) to relate M and Vq after identifying the values of x at iV = 55 numerically. In the regime 
i3 ^ 1 we recover the standard mass scales of inflaton: m/Mpi ~ 6.8 x 10~^ for p = 2 and A ~ 2.0 x 10~^^ for p = 4. 

Since we regard M to be a cutoff scale for the function G{(j),X), the effective theory can be trusted as long as 
H < M. This relation yields the constraint B'^x^ < Mpi/M. For the case B ^ 1 and p = 2, we find that the 
effective theory can be trusted for x < M/m k, (M/10"'^^ GeV)^. For i3 ^ 1 and p = 4, this constraint reduces to 
X < A-i/4(M/Mpi)i/2 « (M/10i4GeV)3/2. 

Finally we note that the scalar propagation speed squared (|175|) . in the regime Sgx 3> Sx, reduces to ~ 2/{3n). 
Also, although the non-Gaussianity parameter /^l"'' is constrained to be small for n = 1, it is possible to have 
I/nl"''| ^ 1 for n ^ 1. It would be of interest to see whether or not such models can be compatible with observations. 



IX. CONCLUSIONS 



In this paper we have studied the observational signatures of chaotic inflationary models with the potential V^cj)) = 
Vo(0/Mpi)P in the context of modiflcd gravitational theories. In Einstein gravity the self-coupling potential ¥{(/)) = 
A(^^/4 is excluded by CMB temperature anisotropy data, while the quadratic potential V{(j)) = m?(tP' /2 is within the 
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2(7 observational contour. Our main aim here has been to clarify how various field couplings present in low-energy 
effective string theory modify the scalar/tensor power spectra generated during inflation. 
We have found a number of new results summarized below. 

• (i) The inclusion of a non-canonical kinetic term uj((l>)X with the exponential coupling — e'^'^/^^p> (/i > 0) 
allows the chaotic inflation models that are in tension with observations to be made compatible with them. 
We have studied the effects of the non-canonical kinetic coupling as well as the nonminimal coupling on the 
inflationary observables and have placed bounds on the strength of the couplings by using recent data from 
WMAP 7yr, BAO, and HST. 

• (ii) In Brans-Dicke theory we have found that the field potential of the form V{(f>) = Vb(0 — MpiY, where p is 
close to 2, can be viable for inflation followed by a successful reheating. We have evaluated the scalar index 
rig and the tensor-to-scalar ratio r for the potential V{(j>) = Vo(0 — Mpi)^ and have shown that r decreases for 
smaller values of the BD parameter wbd ■ The models where wbd is around the order of unity, which includes 
the /(i?) = R + / (6A/^) model, is well within the la observational bound constrained by WMAP 7yr, BAO, 
and HST. 

• (iii) In the presence of the Gauss-Bonnet coupling of the form S,{(t>)G, where ^(0) = ^ge^'^^^^p', we have found 
that the GB coupling with positive /i does not save the self-coupling potential V{(l)) = Xcfi'^/A. For the quadratic 
potential 1^(0) = m^(/i^/2 we have shown that the GB coupling needs to be suppressed (|Jj/es| < 0.1) from the 
CMB likelihood analysis. If iJ. is negative then it is possible to lead to the decrease of both jrig — 1| and r for 
negative ^o, but we require a large coupling constant, such as |^o| ^ 10^", in order to produce a sizable effect on 
the inflationary observables. 

• (iv) In the presence of the Galileon-like self-interaction G(0,X)n<?!), where G{<p,X) = $(<?!))X" and $ cx e^-^/^^p', 
we have expressed the three inflationary observables rig, r, and rit in terms of two slow- roll parameters Sx 
and Sgx- In the regime where the Galileon term dominates over the standard kinetic term {Sgx 3> <5x) we 
have derived analytic formulas for and r in terms of the number of e-foldings N, which recover the results 
obtained for /i = 0. We have shown that, for fi > 0, the Galileon term can lead to the compatibility of the 
chaotic inflationary potentials with current observations. We have confirmed this property for the potential 

= X(j)'^/A by carrying out the CMB likelihood analysis. 

In summary, we have undertaken a unified study of the effects of a number of generalisations to the standard 
inflationary picture as motivated by low-energy effective string theory. We have found that a number of chaotic 
inflationary models which are in tension with observations can be made compatible with them through the addition 
of such terms. The stronger constraints on rig and r expected from the PLANCK satellite will provide an opportunity 
to further test the viability of such scenarios. 
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Appendix A: Fitting function for G-inflation with an exponential coupling 

In this Appendix we present the fitting function for the quartic potential V^(0) = A(/)^/4 in the presence of the 
Galileon- type coupling G ~ — (l/Af^)e^'^/^^p'Ar with /i = 1. Numerically we find the field value giving iV = 55 
before the end of inflation and evaluate Sx and Sqx for several different values of B (B = 10'/^ with i = —32, . . . , 32). 
These slow-roll parameters can be approximated by the following fltting function (found by using the method of least 
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Figure 8: Numerical data points corresponding to the values of Sqx and 5x satisfying the constraint A*' = 55. Each data point 
corresponds to a particular value B = 10'''* with i = —32, . . . , 32. The fitting function used for the CMB likelihood analysis is 
also plotted. 
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We have used this expression in the regime 10^* < Sx < 0.018 for the CMB likelihood analysis in Fig. [T] Finally, in 
Fig. [51 we show both the numerical data and the fitting function 6gx ~ <5gx(<5x)- Since its inverse function, on the 
whole interval, is multivalued, we have used Sx as the independent slow-roll parameter. 
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